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^ . Abstract 

o 

• It is argued that the dominant subleading shape-function contributions to the 

endpoint region of the charged-lepton energy spectrum in i? ^ Xy, I v decays can 
Qh! be related in a model-independent way to an integral over the B Xg'j photon 

spectrum. The square root of the fraction B ^ I v events with charged- 
lepton energy above = 2.2 GeV can be calculated with a residual theoretical 
uncertainty from subleading shape-function effects that is safely below the 10% 
^ ' level. These effects have therefore a minor impact on the determination of iV^fel. 



1. Introduction: One of the most promising strategies for the extraction of the Cabibbo- 
Kobayashi-Maskawa matrix element \Vub\ rehes on the measurement of the inclusive 
semileptonic B X„ I v decay rate in the endpoint region of the charged-lepton energy 
spectrum, which is inaccessible to decays with a charm hadron in the final state . Non- 
perturbative effects can be controlled systematically by using a twist expansion Q and 
soft-collinear factorization theorems @, ^. At leading order in 1/mf,, bound-state effects 
are incorporated by a shape function accounting for the "Fermi motion" of the h quark 
inside the B meson. This function can be determined experimentally from the photon 
energy spectrum in inclusive radiative B — > decays 0. 

Recently, there have been first discussions of the structure of subleading-twist con- 
tributions to the B — > and B I v spectra, which (at tree level) can be param- 
eterized in terms of four subleading shape functions 0. The phenomenological impact 
of these functions on the inclusive determination of has been investigated in [0, |^. 
These authors point out that certain 1/mf, corrections related to chromo-magnetic in- 
teractions appear to be enhanced by large numerical coefficients. They conclude that 
the ignorance about the functional form of the subleading shape functions would lead to 
a significant theoretical uncertainty in the determination of iV^bl, which could only be 
reliably reduced if the lower cut on the lepton energy were taken below the region where 
Fermi-motion effects are important (i.e., below 2 GeV or so). For a value Eq = 2.2 GeV, 
as employed in a recent analysis reported by the CLEO collaboration the resulting 
uncertainty on \ Vub\ was estimated to be at the 15% level 0. While using simple models 
the correction was found to be negative, it was argued that the sign of the effect was 
uncertain in general 

In the present note we explore in more detail the origin of the "enhanced" correc- 
tions found in these papers. Our main point is that the first moments (but not higher 
moments) of the subleading shape functions give a large, non-vanishing contribution 
to the integral over the lepton spectrum even if the lower lepton-energy cut is taken 
out of the endpoint region. This effect corresponds to a calculable correction of order 
Aqqj)/ {mi, AE), where AE = Mb/2 — Eq. The hadronic uncertainty inherent in the 
modeling of subleading shape functions must therefore be estimated with respect to this 
contribution. When this is done, the remaining theoretical uncertainty is found to be 
much less than what has been estimated in [0, |^. We show how the effect of the first 
moments of the subleading shape functions can be isolated and expressed in a model- 
independent way in terms of the photon energy spectrum measured in S — > Xg'y decays. 
We then estimate the numerical effect of the residual higher-twist corrections and find 
their impact on the \Vub\ determination to be small, safely below the level of 10%. 

2. Charged-lepton energy spectrum: The quantity of primary interest to the determi- 
nation of \Vub\ is the normalized fraction of B ^ Xuli' events with charged-lepton 
energy above a threshold Eq chosen so as to kinematically suppress the background from 
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When combined with a prediction for the total B —>■ Xu I v decay rate, knowledge of the 
function Fu{Eq) allows one to turn a measurement of the branching ratio for B X^ I v 
events with Ei > Eq into a determination of \Vub\- 

In the formal limit where the "energy window" AE = Mb/ 2 — Eq is such that 
Aqcd -C AE <C mfc, Fermi-motion effects can be neglected, and the function Fu{Eq) can 
be calculated using the operator product expansion. At tree level the result is 
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where A = Mb — mi,, and 2AE — A = rub — 2Eq is twice the width of the energy window 
in the parton model. The hadronic parameters Ai and A2 measure the fe-quark kinetic 
energy and chromo-magnetic interaction inside the B meson. Note that while the lead- 
ing contribution in (|]) is proportional to the width of the energy window, the power 
corrections are independent of AE. As a result, the relative size of the power correc- 
tions strongly increases as the energy cut Eq is raised toward the kinematic endpoint 
(corresponding to AE — > 0). Although this simple analysis breaks down as AE ~ A, it 
explains that the origin of the large power corrections found in |^, ^ is the kinematic 
suppression of the leading-order term. 

For realistic values of the energy threshold the quantity AE is of order A, and the 
operator product expansion must be replaced by the twist expansion 0, § . At subleading 
order in l/rrib the tree-level expression for Fu{Eq) can then be written as 
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is a combination of the leading and subleading shape functions 0, and the dots denote 
higher-order terms in the expansion. The function ^^,(0;) defined by the second relation 
is related to the normalized photon energy spectrum in B ^ decays, S{E^), by (the 
factor 2 results from the Jacobian duj/dE) 

1 dT{B^ X,7) 
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It is important in this context that the shape of the B — ^ Xg'j photon spectrum is 
largely insensitive to possible effects of New Physics 0, so J^s{^) can be extracted from 
the data in a model-independent way. When we include radiative corrections below, 
S{E^) will still denote the photon energy spectrum, normalized however on an interval 
^min < < Mb/2, with i?™™ sufficiently small to be out of the shape-function region. 

The combination of subleading shape functions remaining in the last line of pa- 
rameterizes chromo-magnetic interactions in the B meson. The moment expansion of 
these functions yields 
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where p2 is a 5-meson matrix element of a local dimension-6 operator. In the limit 
where AE ^ A, only the first moment yields a non-zero contribution to the function 
Fu[Eq) in (0), because the weight function under the integral is linear in u (to first order 
in 1/mb). On the other hand, near the endpoint of the lepton spectrum all moments 
of the shape functions become equally important 0, ^ . In between these two extremes 
there is a transition region, where only the first few moments of the shape functions 
give significant contributions. Theoretical studies of the photon spectrum in i? — > Xgj 
decays have shown that this transition region corresponds to values Eq ~ 2.0-2.3 GeV 
(for yet lower values, Fermi- motion effects become unimportant) To account for the 
effect of the first moment we define a new subleading shape function 
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whose normalization and first moment vanish, and whose contribution to the quantity 
Fu{Eq) therefore vanishes for AE ^ A. Inserting this definition into relation (^, and 
using that ^-"5(0;) = /(cu) + . . . to leading order in 1/mf,, we obtain from (□) 
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Taking into account the known 0{as) corrections to the leading term in the twist expan- 
sion fiy, |ll|, and rewriting the contribution involving J-'s{uj) as a weighted integral over 
the normalized photon energy spectrum in 5 — > Xsj decays, we get our final result^ 
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^Using an integration by parts, this result can be rewritten as a weighted integral over the fraction 
Fs(E) oi B ^ Xsj events with photon energy above E, normalized such that i^s(£'™'") = 1- 
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and the subleading shape-function contribution 
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The factor 2 in front of Asl(-E'o) in ® is inserted so that AsL{E())/mf, is the subleading 
shape-function correction to \Vub\- We stress that, by definition, Asl(-Eo) is a parameter 
of order Aqcd that vanishes for AE ^ A. It is thus a true measure of shape-function 
effects. On the contrary, the power corrections studied in ^] arise predominantly from 
the X2/1TLI correction to the weight function. 

The expression for the perturbative coefficient fcpert in ( P^ can be obtained from the 



results of [Bl 121. It reads 



wi^r) = -f -^1-'^-^!-^+ E ^j^^M^)^ (12) 

where 5=1 — -E™"/ {E^) depends on the lower boundary of the energy interval used to 
normalize the B Xs'~f photon spectrum, Ci{fi) are leading-order Wilson coefficients in 
the effective weak Hamiltonian for B Xg'y transitions, and the functions fij{S) can be 
found in [01. In the definition of 6 we use the central value of the CLEO result for the 



average photon energy above 2GeV, {E^) = (2.346 ± 0.034) GeV |T^, as a substitute for 
mb/2. 

3. Numerical results: The value of the coefficient fcpert is sensitive to the choice of 
the renormalization scale /i and the value of the quark-mass ratio mc/rrih used in the 
evaluation of charm-quark loops. We take fib = rni,{mh) = 4.2 GeV as our central value 
for the renormalization scale and vary /z between /if,/2 and 2/ife. We use a running charm- 
quark mass to evaluate the loop functions [|14|, taking mc{fJ')/mf,{fi) = 0.23 ± 0.03. The 



results for kpert corresponding to two different choices of ii^™" are 

^ ^ 0.07±0.10; E--=1.5GeV, 

^ -0.20 ± 0.08 ; Ef"" = 1.75 GeV . ^ ^ 

Since the effect of this correction is very small, one should not consider the small vari- 
ation of fcpert as a measure of the perturbative uncertainty in the weight function ([10|). 
Typically, we expect O(a^) corrections to contribute at the level of 5% of the tree-level 
term. A corresponding uncertainty will be included in our numerical analysis below. 
The power correction to the weight function in ( pUj ) may be rewritten as 

8A2 m'i,, — m%, , , ,^ 

^ ^ 0.044 . (14) 
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Alternatively, using A2 = (0.12 ± 0.02) GeV^ and mt = (4.72 ± 0.06) GeV we obtain 
the value 0.043 ± 0.007, which will be used in our numerical analysis. The size of this 
correction is not anomalously large; however, its impact is significant because it competes 
with terms proportional to the small difference (1 — Eq/E^). In the endpoint region this 
difference scales like AqcD/'^ft, and so the X2/^l term is of relative order l/rrib. 

Our final focus is on the subleading shape-function contribution Asl{Eo) defined in 
(O). Little is known about the function s{u)) except that its normalization and first 

(s) 

moment vanish, and that its second moment, = —p2, is given by a hadronic matrix 
element expected to be of order (0.5 GeV)'^ with undetermined sign. As a result, the 
functional form and sign of Asl(-Eo) cannot be predicted at present. However, the fact 
that Asl(-Eo) must approach zero as Eq is lowered to a value of about 2 GeV (below 
which shape-function effects from higher moments are irrelevant) ensures that its impact 
on the determination of \Vub\ is small. To substantiate this claim we investigate several 
models for the subleading shape function in more detail. For the leading-order function 
we take the ansatz |^ 



fi^) = j^9a{x) , with X 
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1 + T > 
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where ga{x) = [a'^/T{a)] e~"^. The parameter a must be larger than 1 and is fixed 
so that the second moment of /(cu) equals — Ai/3 |@|, yielding a = — SA^/Ai. We assume 
that the subleading function s{uj) is finite everywhere in the interval —A < u < 00, but 
we do not require that this function vanish at the endpoint. 

The model functions adopted in ||^ are such that s{u) is set to zero, and so Asl(-Eo) 
vanishes by construction. The model functions used in [H correspond to the ansatz 
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where the lower bound on the parameter b is enforced by the requirements that s{uj) be 
finite at the endpoint u = —A and have vanishing normalization and first moment. A 
property of this model is that also the second moment of s{uj) vanishes. Three alternative 
choices for s{uj) with non-zero second moment M2'^^ are 
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Figure |1| shows results for Asl(-Eo) obtained in the various models, using A = 0.5 GeV, 
Ai = —0.3 GeV^, and Mg*"* = (0.5 GeV)^ as input parameters, and varying the parameter 
b over a wide range of values. Although the details of the subleading shape function s{uj) 
are rather different in the four cases, all models exhibit the same general features. While 
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Figure 1: Model predictions for the subleading shape-function correction Asl(-Eo) 
as a function of the cut Eq. For each model, the parameter b is varied between the 
minimal allowed value (red) and 10 (blue) in steps of 1. The sign of Asl(-Eo) is 
undetermined. 



Asl(-Eo) can be large close to the kinematic endpoint, it takes values of order A for 
£^0 ~ 2.35 GeV and quickly decreases as Eq is lowered below 2.3 GeV. For Eq = 2.2 GeV 
we find values of Asl(-Eo) of at most 130 MeV (model 2), corresponding to a power 
correction to the extraction of \Vub\ of less than 3%. Although our choice of model 
functions is meant as an illustration only, we believe the rapid decrease of Asl(-E'o) 
for Eq < 2.3 GeV is a general result. It appears to be extremely unlikely that with 
a reasonable shape of s{uj) and a natural size of the second moment M2 the power 
correction AsL(-Eo)/^fe could be as large as 10%. 

4- Conclusion: In summary, we have studied the impact of subleading shape functions 
on the determination of \Vub\ from the combination of weighted integrals over energy 
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Table 1: Illustrative theoretical predictions for the fraction Fu{Eq) oi B ^ Xu I v events 
with charged- lepton energy Ei > Eq, assuming a perfect measurement of the B — > Xgj 
photon spectrum (see text for explanation). 



Eo [GeV] 


LO 


NLO 


l/nih 


total 


residual error 


2.0 


0.271 


0.041 ±0.014 


-0.040 ±0.006 


0.273 ±0.015 


±0.003 


2.1 


0.195 


0.033 ±0.010 


-0.037 ±0.006 


0.191 ±0.011 


±0.005 


2.2 


0.126 


0.024 ±0.006 


-0.033 ±0.005 


0.117 ±0.008 


±0.006 


2.3 


0.068 


0.015 ±0.004 


-0.026 ±0.004 


0.057 ±0.006 


±0.008 



spectra in inclusive B I v and B — Xp( decays. We have argued that for a 

lower energy cut Eq = 2.2 GeV as employed in a recent CLEO analysis one is in a 
transition region, where Fermi-motion effects are dominated by the first few moments 
of the leading and subleading shape functions. The dominant power correction (the 
only one that remains when the cut is lowered below about 2 GeV) results from the 
first moment of the subleading shape function, which is known in terms of the hadronic 
parameter A2. 

Our main result is given in (^ and (0). To exhibit its features, let us assume that a 
perfect measurement of the B —>■ Xs'~f photon spectrum is available in the energy range 
above i?™™ = 1.5 GeV. (For the purpose of illustration, we use a fit to the CLEO data 
in ||1^.) We then calculate the fraction of i? — > X^ I v events with charged-lepton energy 
above Eq for different values of the cut. The results are summarized in Table |I]. Columns 
2, 3, and 4 show the contributions from the tree-level term, the 0{oi.^ corrections, and the 
power correction to the weight function in ([T0|), including theoretical uncertainties from 
input parameter variations as detailed above. The next column shows the total result, 
while the final column gives an estimate of the residual uncertainty from subleading 
shape-function effects, as parameterized by the term 2AsL(-Eo)/mfe in (H). We show 
the largest uncertainty obtained in the four classes of models considered earlier. We 
observe that the power correction to the weight function has a significant impact, which 
as anticipated is by far the dominant effect of subleading shape functions. For Eq = 
2.2 GeV, the power correction leads to a reduction of the predicted value for Fu{Eq) by 
(26 ± 6)%, corresponding to a 13% enhancement of the extracted value of \Vub\- This is 
in good agreement with the estimate given in [0. 

The most important implication of our analysis is that subleading shape-function 
effects do not entail a significant limitation on the extraction of \Vuh\- This assessment 
differs from the conclusion reached in §| , where is was argued that these effects could 
not be controlled reliably unless the cut Eq could be lowered outside the shape-function 
region. The new element of our analysis is that we identify the first moment of the 
subleading shape-function as the dominant source of power corrections and show how 
its contribution can be expressed in terms of an integral over the B photon 



7 



spectrum. We have estimated the residual uncertainty on \Vub\ from subleading shape- 
function effects by using four different classes of model functions and found corrections 
of at most 3% (with Eo = 2.2 GeV). The smallness of this effect can be understood on 
the basis that it is a power correction of the form Ash{Eo)/mb with a hadronic parameter 
Asl(-Eo) = O(Aqcd) that vanishes as Eq is lowered below about 2 GeV. We thus conclude 
that, very conservatively, the residual uncertainty on \Vub\ is less than 10%. 

The main result of this letter is the new expression for the weight function in ([T0|), 
which now includes the leading power correction. Perhaps the largest uncertainty in this 
method for determining \Vub\ is due to (largely unknown) corrections from violations of 
quark-hadron duality, and from spectator-dependent effects such as weak annihilation 
and Pauli interference [^, [T^ (see also where a 6-8% correction on \ Vub\ was obtained 
for Eq = 2.2 GeV using a simple model for spectator effects). In these references, several 
strategies have been developed that could help to determine the magnitude of these 
corrections using experimental data. 
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